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Abstract 

jL We present the results of an optical and X-ray monitoring campaign on the short-period massive SB2 binary HD 152218. Combining 
< l our HiRes spectroscopic data with previous observations, we unveil the contradictions between the published orbital solutions. In 

particular, we solve the aliasing on the period and derive a value close to 5.604 d. Our eccentricity e = 0.259 ±0.006 is slightly lower 

than previously admitted. We show that HD 152218 is probably undergoing a relatively rapid apsidal motion of about 3° yr^ 1 
1 'and we confirm the 09IV+09.7V classification. We derive minimal masses of 15.82 ± 0.26 Mq and 12.00 ± 0.19 M© and constrain 

rhe radius of the components to Ri = 10.3 ± 1.3 R© and R2 = 7.8 ± 1.7 R©. We also report the results of a yLMM-Newton 
^ monitoring of the HD 152218 X-ray emission throughout its orbital motion. The averaged X-ray spectrum is relatively soft and it 

is well reproduced by a 2-T optically thin thermal plasma model with component temperatures about 0.3 and 0.7 keV. The system 



in 



presents an increase of its X-ray flux by about 30% near apastron compared to periastron, which is interpreted as the signature of 



£-««. an ongoing wind-wind interaction process occurring within the wind acceleration region. 
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.h 1' Introduction 

X: 

^ . Early-type stars of spectral type O are character- 
" ized by strong stellar winds. Though not as extreme 
as those of their evolved descendants, the Wolf-Rayet 
stars, these winds combine terminal velocities of a few 
thousand kins -1 and important mass-loss rates (about 
10~ 7 - 1(T 5 M Q yr" 1 ) that significantly affect both the 
surroundings of the star and its evolution. In a binary 
system, it is expected that the winds from the two stars 
collide, producing a density enhanced region known as the 
wind interaction zone. The shock-heated plasma within 



* Based on observations collected at the European Southern Obser- 
vatory (La Silla, Chile), at the Cerro Tololo Inter-American Observa- 
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Member States and the USA (NASA). 
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**Visiting Astronomer, CTIO, National Optical Astronomy Obser- 
vatories (NOAO). NOAO is operated by the Association of Univer- 
sities for Research in Astronomy, Inc. under contract with the Na- 
tional Science Foundation. 



this zone is expected to produce an additional contribution 
to the X-ray emission from the early-type system and, in- 
deed, the early-type binaries are known to be statistically 
X-ray over lum inous compared to single stars of the same 
spectral type ( Chlebowski fc Garmanv . 199ll ). This extra 
X-ray emission might display phase-locked modulations, ei- 
ther due to a change of the absorption properties along the 
line of sight, or to a modulation of the shock strength due 
e.g. to a varying separation between the components in an 
eccentric binary. The properties of the shock thus strongly 
depend both on the geometry of the binary system and on 
the characteristics of the individual winds. However these 
properties often cruelly lack accurate constraints. In this 
framework, our team is involved in a long-standing effort 
to study early-type stars in a number of open clusters. 

The present work focuses on the short-period binary 
HD152218, located in NGC 6231 at the core of the ScoOB 1 
assoc iation. This star is a long known SB2 system dStruve , 
19441) Previous papers (iHill et al.l.ll974tlLevato fc Morrel , 



19831: IStickland etail Il997t iGarria fc Mermilliodl . 12001 ) 



however reported conflicting results and we thus decided to 
re-appraise the orbital properties of this apparently well- 



Table 1 

Journal of the spectroscopic observations of HD 152218. Column 1 lists the Heliocentric Julian Date (in format HJD— 2 450 000) at mid- 
exposure. The next three columns present the phases (calculated from the Hel-line SB2 orbital solution of Table |4]l and the averaged primary 
and secondary radial velocities computed in the respective systemic velocity frames (see Sect. I3.3|l . Instrumental setups are given at the 
bottom of the table. 



HJD 

(days) 


0HeI 






HJD 0Hd 

(days) 






RV\,l - 7a, l RV\,2 - 7a, 2 
(bus- 1 ) (kms" 1 ) 


RV\,i - 7a, l RV\,2 - 7a, 2 
(kms -1 ) (kms -1 ) 


621.673° 


0.406 


-25.00 


-29.31 


1302.901 c 0.968 


-2.09 


0.46 


622.642° 


0.579 


63.69 


-133.73 


1304.902 c 0.325 


-82.39 


134.85 


622.755° 


0.599 


80.30 


-123.20 


1323. 853 c 0.707 


132.51 


-175.86 


623.648° 


0.759 


151.20 


-200.79 


1327.796 c 0.411 


-11.80 


-9.25 


623.762° 0.779 


163.30 


-188.40 


1328.905 6 0.609 


82.27 


-124.90 


624.590° 


0.927 


74.35 


-116.83 


1329.855 6 0.778 


147.17 


-199.78 


624.722° 0.950 


13.36 


9.05 


1331.882° 0.140 


-182.39 


213.02 


625.626° 


0.112 


-158.93 


209.21 


1669.812 c 0.442 


-6.42 


-3.87 


625.734° 


0.131 


-171.70 


233.10 


1670.804 c 0.619 


96.65 


-114.87 


626.591° 


0.284 


-104.57 


166.73 


1670. 918 c 0.639 


108.37 


-132.91 


626.716° 


0.306 


-85.00 


141.40 


1671. 828 c 0.802 


153.10 


-200.44 


1299.883 c 


0.430 


-7.56 


-5.01 


1672. 865 c 0.987 


-12.61 


-10.06 


1300.882 c 


0.608 


94.15 


-107.22 


1673. 878 c 0.168 


-168.93 


211.74 


1301. 883 c 


0.787 


152.44 


-201.56 


2037.812 c 0.110 


-172.15 


221.33 



studied system. The aim is to bring its orbital and physical 
parameters on firm ground. As mentioned above, accurate 
ephemeris and a detailed knowledge of the orbital and 
physical properties of a system are indeed crucial ingredi- 
ents for interpreting the X-ray data and for uncovering the 
possible signature of a wind interaction. 

The remaining of this paper is organised as follows. The 
next section describes the observational material and the 
data reduction processes. Sect.[3]aims at unveiling the dis- 
crepancies of the previously published orbital solutions. 
Sect. |4] discusses the physical properties and evolutionary 
status of HD 152218 and, in Sect. El we analyse the X-ray 
observations in the light of the upgraded orbital solution. 
Finally Sect. [5] provides a summary of the present work. 

2. Observations and data reduction 
2.1. Optical spectroscopy 



a Loral 2688 x 512 pixel CCD with a pixel size of 15 x 15 ^m. 
The slit width was chosen to achieve a nominal resolving 
power of 70000-80000. The effective resolving power as de- 
rived from the FWHM of the lines of the ThAr calibration 
exposures is 65000-75000. Typical exposure times were of 
the order of 30 minutes and the average signal-to-noise ra- 
tio (SNR) is about 140. The observed wavelength domain 
is centered on the He 1 14471 or He n 14686 line and is ~ 45 A 
wide. 

Another set of three echelle spectra over the range 
3850 to 5790 A was obtained with the Bench-Mounted 
Echelle Spectrograph (BME) attached to the 1.5 m CTIO 
Ritchey-Chretien Telescope, during a 5-night run in June 
1999. Forty nine orders were observed using the KPGL2 
316 lines mm -1 grating as a cross-disperser. The detector 
was a Tek 2048 CCD with 24 /an pixels. The slit width was 
set to 70 fim corresponding to a resolving power of 45000. 
Exposure times ranged from 45 to 60 min, yielding a SNR 
in the continuum between 40 and 50. 



The present work is based on 42 high-resolution spec- 
tra of HD 152218 obtained at the European Southern Ob- 
servatory (ESO, La Silla, Chile) and at the Cerro Tololo 
Inter- American Observatory (CTIO). These were acquired 
during 35 different nights distributed over 8 runs between 
May 1997 and May 2004. The journal of the observations 
is presented in Table [TJ 

In May 1997, five (resp. six) high resolution spectra of the 
He 1 14471 (resp. He n 14686) line were obtained with ESO's 
1.4 m Coude Auxiliary Telescope (CAT) at La Silla, using 
the Coude Echelle Spectrometer (CES) equipped with the 
Long Camera (LC). The detector used was ESO CCD#38, 



Between April 1999 and May 2002, we collected 20 
echelle spectra covering the whole optical range (~3750- 
9200A) using the Fiber-fed Extended Range Optical Spec- 
trograph (FEROS), an echelle spectrograph mounted at 
the ESO 1.5m telescope at La Silla. In May 2003 and May 
2004, respectively three and five other FEROS spectra 
were obtained at the ESO 2.2m telescope, at La Silla too. 
The detector was a 2k x 4k EEV CCD with a pixel size of 
15/im x 15/im. The spectral resolving power of FEROS 
is 48 000. Typical exposure times ranged, at the ESO 1.5m 
telescope, from 10 to 20 min according to the weather 
conditions, resulting in typical SNR around 150. At the 



Table 1. Continue 
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83 
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34 
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10 
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a CES + ESO-CAT; b BME + CTIO 1.5m ; c FEROS + 
ESO 1.5m; d FEROS + ESO/MPG 2.2m 

ES0 2.2m telescope, the SNR is about 200 in May 2003 
and about 250 in May 2004. Typical exposure times at the 
ESO 2.2m were of 13 min. 

The CES data were reduced in a standard way using 
the MIDAS package supported by ESO. The spectra were 
rectified by means of an instrumental response curve built 
from the observations, obtained under similar conditions, 
of a metal-poor 'reference' star (HD 203608: F6V). Finally 
the spectra were normalized by fitting a low order polyno- 
mial to the continuum. 

The BME data were reduced using the IRAH3 pack- 
age and following the recommendations of the BME User's 
Manual. The pixel to pixel variations were removed using 
flat field exposures taken with a very bright light source and 
a diffusing screen placed inside the spectrograph (so-called 
milky flats). A first rectification of the extracted echelle or- 
ders was carried out with the projector flat exposures. The 
spectra were then normalized by fitting a low-order poly- 
nomial to the continuum. 

FEROS data were reduced using an improved version 
of the FEROS context (see details in Sana et all . [2003) 
working under the MIDAS environment. We mainly used 
the normalized individual orders. However, when doubts 
about the quality of the normalisation arise, we checked 
our results with the automatically merged spectrum. 

2.2. X-ray observations 



Newton European observatory (jJansen et all 120011 ). The 
campaig n, of a total duration of about 180 ks, is described in 
length in I Sana et al.l (|2006a ) together with the data reduc- 
tion and source identification. We only give a brief overview 
here, focusing on the additional elements specific to HD 
152218. 

The campaign actually consisted of six separate obser- 
vations spread over 5 days. All three E PIC instruments 



(jStriider et all l2001t iTurner et all |2001() were operated in 



the Full Frame mode together with the Thick Filter to 
reject UV/optical light. The field of view (FOV) of the 
EPIC instruments was centered on the colliding wind bi- 
nary HD 152248, from which HD 152218 is separated by 
about 6.'8. Due to the brightness of the cluste r obje cts in 
the FOV, the Optical Monitor (jMason et al.l . l200lh was 
switched off throughout the campaign. Fig. [1] presents a 
view of the combined MOS images around HD 152218 while 
Table [5] provides the journal of the X-ray observations of 
the system. 

HD 152218 unfortunately fell on a gap of the pn instru- 
ment and, as in lSana et al.1 ( 2006b ). we only considered the 
EPIC MOS data. To obtain the averaged count-rates dur- 
ing the different pointings , we adopte d the so urce positions 
given in the catalogue of ISana et al. I (<2006bl) and we per- 
formed a psf (point spread function) fit using the SAS task 
emldetect. As in the previous works, we restrained our sci- 
entific analysis to the 0.5-10.0 keV band and we adopted 
three energy sub-ranges: asoft (Sx) [0.5:1.0 keV], amedium 
(M x ) [1.0:2.5 keV] and a hard (H x ) [2.5-10.0 keV] band. 

Located slightly offset from the cluster core, HD 152218 
is reasonably well isolated in the X-ray images. To extract 
the source spectra, we used an extraction radius of 20'.'2, 
centered on the source position for all six pointings and 
for the two EPIC MOS instruments. The different back- 
grounds were estimated from source-free regions located 
on the same CCD detector as HD 152218 (CCD #3 and 5 
respectively for the MOS1 and MOS2 instruments). The 
adopted source extraction region is shown in Fig. [1] Dis- 
plays of the cor responding backgr ound extraction regions 
are presented in lSana et al.l (|2006d l together with their po- 
sitions and sizes. 

Using these regions, we extracted the source spectra from 
the two MOS instruments and for the six pointings individ- 
ually. We also extracted the merged spectra in each instru- 
ment, thus combining the six observations of HD 152218. 
For this purpose, we built the corresponding Ancillary Re- 
sponse Files (the so-called arf files) using the SAS task 
arfgen. We adopted the Redistribution Matrix Files (rmf) 
provided by the SOC. The spectra were finally binned to 
reach at least 10 counts per bin. The X-ray data will be 
discussed in Sect. [5] 



In September 2001, the NGC 6231 open cluster was the 
target of an X-ray observing campaign with the XMM- 

1 IRAF is distributed by the National Optical Astronomy Observa- 
tories. 
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Table 2 

Journal of the XMM- Newton observations of HD 152218. The Julian Date (JD) at mid-exposure is given in Col. 2. Cols. 3 and 4 list the 
effective exposure times for the EPIC MOS1 and MOS2 instruments while Cols. 5 and 6 report the background-subtracted, vignetting and 
exposure-corrected count-rates in the different instruments. The last column provides the orbital phase of HD 152218 for each XMM- Newton 
observation at mid-exposure, according to the Hel-line ephemerides given in Table [4] The phase extension of the pointings is about 0.06. 
Obs. 2 and 4 were however affected by high background events that reduced the nominal exposure time by about one third. 



Obs. # JD Effective duration (ks) Count rates (10 3 cnts 1 ) Phase </>Hel of 
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26.0 


24.3 


12.5 ±0.9 


11.6 ±0.9 
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5 


2161.774 


30.9 


31.0 
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13.4 ±0.9 


0.226 


6 
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32.9 


32.8 
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Fig. 1. Combined EPIC MOS1 + MOS2 image in the vicinity of HD 
152218. The adopted source extraction region is shown. Neighbouring 
X-ray s ources are labelled using the internal numbering of lSana et al.l 
(20063). 



(FWHM) of the fitted Gaussian profiles. We thus computed 
the mean FWHMs of the primary and secondary lines ac- 
cording to the measurements performed near quadrature. 
Considering the properties of the studied lines, we then de- 
cided either to fix the line FWHMs to their mean value or 
to request the primary and secondary lines to have equal 
FWHMs. The resulting radial velocities (RVs) were as- 
signed a lower weight in the determination of the orbital so- 
lution (see Sect. 13. 3| ). In order to get a reliable solution, we 
considered if different lines: He 1 114026, 4388, 447i, 4922, 
5016, 5876, 7065, He nil 4542, 4686, 5412 and O m 15592. 
Note that the HeiiAA4542 and 5412 lines were rarely de- 
blended and show a faint secondary component only. For 
these lines, we thus focused on the primary orbital motion. 
We also note that the O in 15592 line is only seen in the 
primary spectrum. Finally, to compute the RVs associated 
with the measured Doppler shif ts, we adopted the effec- 
tive wavelengths for O-stars from Conti et al.l (jl977l ) below 



4800 A and from lUnderhiill (|l994h above 



3. Unveiling the correct orbital solution 

3.1. Optical spectrum 

The spectrum of HD 152218 is clearly dominated by the 
usual Balmer, Hei and Hen lines seen in absorption (Fig. 
[2]). Beside these, the spectrum reveals metallic lines of N in, 
Si in, Si iv, Civ and Oni that are typical of O-type stars. 
Only the C in 15696 line is seen in emission. Most of the 
observed stellar lines are double (see e.g. Fig.[3j) and display 
clear Doppler shifts. The HD 152218 spectrum also reveals 
a number of diffuse interstellar bands (DIBs), as well as 
interstellar absorptions due to Nai, Cai, Can, CH and 
CH+. 

We measured the Doppler shifts by fitting a sum of Gaus- 
sians to the different line profiles. Whenever the compo- 
nents get closer to conjunction, the disentangling of their 
spectral signature was more difficult. To get reliable re- 
sults we had to constrain the full width at half maximum 



3.2. Period determination 



Our period search relies on the method of lLafier &: Kinman 
( 19651 LfeK) an d on the Fourier-an alysis technique of 



Hec k et al.l (119851 IHMM85I . see also iGosset et al.l 2001 



for comments) applied on the different data sets of Ta- 
ble [3J Using the L&K algorithm, the obtained periods 
are in excellent agreement, yielding a mean value of 
5.60380 ± 0.00005 d. The results from the Fourier anal- 
ysis are more conflicting. The Hell 4471 and He ill 4686 
data sets yield a period of 5.604 days while the other lines 
indicate 5.519 days. An inspection of the corresponding 
periodograms revealed that the period determination is 
strongly affected by a 1-year aliasing, leading to a possi- 
ble confusion between orbital frequencies of 0.1784 and 
0.1812 d _1 . These latter values correspond respectively to 
the two period values quoted here above. From Fig. 0] (left), 
it is clear that the period can not be reliably constrained 
from the sole FEROS data. The inclusion of the CES and, 
to a lesser extent, of the BME data clearly improves the 



.1 




Fig. 2. Spectrum of HD 152218 as observed with the FEROS spectrograph on HJD 2 453 134.723 (<p ~ 0.85). At this orbital phase, the 
spectrum of the most massive star is shifted towards the red. The most prominent lines arc indicated at their rest wavelengths. 
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Fig. 3. He 1 14922 (left) and He 1 15875 (right) lines at selected phases. 
The spectra were shifted along the vertical axis according to the value 
of their phase (right-hand scale). The secondary spectral signature 
is clearly identified on most spectra but the blended ones. 
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Fig. 4. Left panel: HMM85 periodo gram (i.e. amplitude spec- 
trum) associated with the FEROS He 1 14471 primary data set. Peri- 
odograms associated with other lines are very similar. Right panel: 
idem, using the CES and BME measurements in addition to the 
FEROS data. This latter data set favors a period of 5.604 d, corre- 
sponding to v K, 0.1784 d _1 . 

situation (Fig. [4] right) but still does not allow to draw a 
firm conclusion. 

Our data sets consist of spectra spread over 2500 days 
for the He 1 14471 and He II 14686 lines or over 1800 days for 
the other lines. These observational time bases correspond 
respectively to natural peak widths Av of 0.0004 d _1 and 
0.0006 d — 1 . The full width at half maximum of the ac- 
tual peaks in the periodogram is in agreement with these 
theoretical values. Adopting an uncertainty of one tenth 
of the peak width, we respectively obtained ap = 0.0013 
or 0.0017 d. These values however only account for the 
uncertainty due to the finite width of the peak, but do not 
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Fig. 5. Left panel: IHMM85I periodo gram corresponding to the 
IUE+optical data sets used by ISLP97I . Right panels: idem, com- 
bining our data with published data sets. 

consider the possible confusion between different aliases. 
To sol ve the aliasing p r oblem, w e tur ned to the litera- 



ture. As Kill et all (|1974 IHCB74D and iLevato fc Morrell 



derived quite different period values 



pared to the ones quoted above, we first applied our pe- 
riod search methods using their data sets and found a very 
large number of aliases covering a frequency range corre- 
sponding to period values between about 5 and 6.7 days. 
Their proposed period values were thus much more poorly 
constrained than what their error - bars su ggested. On the 
other hand, IStickland et al.l (|l997l ISLP97I ) obtained a pe- 
riod P w 5.604 d, ver y similar to one of our aliases. We 
used both the L&K and HMM85 algorithms on the different 
data sets considered bv lSLP97l and we re-derived very sim- 
ilar values, except for the related uncertainties. Adopting a 
typical error corresponding to one tenth of the periodogram 
peak width, we obtained values about 6 times larger. The 
periodogram (Fig. [5j left panel) also indicates a possible 
confusion with a neighbouring alias at v w 0.183 d . The 
period value is again not so clearly established. We how- 
ever note that this latter alias is not present in our data set. 
By combining the different data sets, we could thus hope 
to definitely constrain the period. 

To combine our data with other observational data sets, 
we first averaged the RV measurements obtained from the 
diffe rent lines quoted in Tabled. Then, applying the L&K 
and IHMM85I period searches to all the available RV mea- 
surements, we derived values of P = 5.6039 and 5.6040 d 
respectively. The corresponding Fourier periodogram is dis- 
played in Fig. [5] (right) and clearly shows that the orbital 
peak is now well isolated. The time span of the whole ob- 
servational set is ^21 915 d, yielding a theoretical peri- 
odogram peak width of Av = 4.6 x 10~ 5 d _1 . The observed 
width of the peak agrees with this value and yields AP = 
1.4 x 10~ 3 d. We adopted a final uncertainty corresponding 
to one tenth of the peak width, i.e. ap = 1.4 x 10~ 4 d. 



2 In this first approach, we did not correct the individual line RVs 
for their systemic velocities but we rather performed a direct average. 
Indeed the systemic velocities should be obtained from RV curve 
fitting. They thus require an first accurate estimate of the period 



Table 3. Orbital solutions for HD 152218 as derived from different lines. The usual notations have been adopted. s y /s x gives the adopted secondary to primary ratio of the uncertainties 
on the measured RVs (see Sect. 13. 31 . The quoted error-bars are the 1-cr uncertainties. 
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He 1 15876 5.60410 


1.6 


0.271 ± 0.006 104.5 ± 1.5 162.4 ± 1.1 221.4 ± 1.5 


-21.9 ±1.2 


-21.1 ± 1.5 


4.8 


He 1 17065 5.60410 


1.7 


0.255 ± 0.008 101.0 ± 2.2 165.1 ± 1.6 214.0 ± 2.1 


-17.1 ± 1.7 


-11.7±2.0 


6.4 


Heiil4542 5.60400 




0.288 ± 0.014 101.3 ± 2.2 160.9 ± 2.1 


-21.1 ± 1.2 




6.5 


He ii 15412 5.60419 




0.302 ± 0.011 104.0 ± 1.6 161.4 ± 1.7 


-16. 5± 1.0 




5.2 


in 15592 5.60409 




0.274 ± 0.012 106.1 ± 2.0 162.7 ± 1.9 


-17.1 ± 1.7 




6.0 



3.3. Orbital elements 

Now that the period is reliably constrained, we can de- 
rive the orbital parameters of the system by fitting the ob- 
served RV-curves assuming a Keplerian motion for the two 
components. We first computed a series of orbital solutions 
using the RV sets associated with the different absorption 
li nes listed in TableBl For SB1 lines, we used the algorithm 



of I Wolfe et alj (|19671 ) in its original form and assigned the 



same weight to each measurement. For the data sets associ- 
ated with SB2 lines, we used a modified version of the algo- 



rithm, adapted to SB2 systems as explained in lSana et al 
(|2006al ). and we only considered those data points for which 
the two components could be measured with a reasonable 
accuracy. We assigned a half weight to those RV points for 
which the FWHMs had to be constrained (see Sect. 13. 
Because of the lower SNR of the BME data, we gave them 
a weight of one fourth while computing the orbital solu- 
tion. In the observed wavelength domain, the BME mea- 
surements were further not included for the fainter lines 
nor for the lines more difficult to disentangle (Hei AA4388, 
5016, He II 14542, Oin 15592). In our derived solutions, we 
adopted a period value and a relative primary to secondary 
weight ratio (s y /s x ) that yield the lowest x 2 . These best 
solutions, including a reappraisal of the periocTH. are pre- 
sented in Table [3] together with the corresponding root- 
mean-square (r.m.s) residuals of the fit. 

In a next step, we computed the mean RVs of the SB2 
He I- lines quoted in Table [3] as well as of the sole primary 
lines. For this purpose, we first shifted the individual RVs 
to a common frame, by subtracting the individual systemic 
velocities. The resulting mean RVs are listed in Table [TJ 
The orbital solutions computed using the averaged primary 
(SB1) and the averaged Hei-line (SB2) RV measurements 
are given in Tabled] together with the physical parameters 
of the system. Fig. [5] displays the RV curves correspond- 
ing to the averaged SB2 solution. The latter solution will 
be adopted throughout the rest of the paper to constrain 
the physical properties of the HD 152218 components (see 
Sect. SI). Corresponding primary and secondary apparent 
systemic velocities, as obtained from a weighted mean of 
the values quoted in Table 02 are tT = —18.7 ± 0.6 kms -1 
and yi = -19.2 ± 0.7 kms" 1 . 

Our final SB2 solution is si gnificantly different from the 
most r ecent determination by lGarcia fc Mermilliodl ( 2001 , 
GMOll ) who obtained P ~ 4.9 d and e ~ 0.4. We how- 



ever mostly re cover, w ithin the error-bars, the results of an 
older study bv lSLP97 Indee d, our new values fall within 3- 
a from the results of lSLP97l . though with a slightly smaller 
eccentricity and larger RV-curve semi-amplitudes. Only the 
longitude of t he peri astron seems to differ significantly: 
80?6 ± 2?1 for [SLP97J versus 104?2 ± 1?6 for our solution. 
This suggests the presence of an apsidal motion of the or- 
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Fig. 6. HD 152218 RV-curves corresponding to the He I solution of 
Table [4] Note that the respective systemic velocities have been sub- 
tracted. RV measurements listed in Table [T] have been overplotted. 
Filled and open symbols are for the primary and secondary compo- 
nents respectively. Different symbols indicate different instruments 
used: CES, squares; BME, triangles; FEROS, circles. 

Table 4 

Orbital and physical parameters of HD 152218 as deduced from 
two different data sets: the averaged RVs computed over selected 
primary SB1 lines (Col. 2) and over He I SB2 lines (Col. 3). The usual 
notations have been adopted. To is the time of periastron passage 
and is adopted as phase zero (i.e. <f> = 0.0). Note that the primary 
and the He I lines solution were computed in the systemic velocity 
frame. 





Prim. 


He I lines 


p (d) 


5.60404 


5.60391 


s y / s x 


n. 


1.5 


mi/ 7712 


n. 


1.319 ± 0.014 


e 


0.275 ± 0.008 0.259 ± 0.006 


u> (°) 


104.2 ± 1.6 


104.0 ± 1.6 


Tg (HJD 


3196.493 


3197.229 


-2 450 000) 


±0.024 


± 0.025 


K ± (kmr 1 ) 


162.4 ± 1.1 


162.2 ± 1.2 


K 2 (kins" 1 ) 


n. 


213.9 ± 1.5 


71 (kms -1 ) 


0.7 ± 0.9 


0.9 ± 1.2 


72 (kms -1 ) 


n. 


-0.6 ± 1.4 


ai sin i (Rq) 


17.28 ± 0.13 


17.34 ± 0.13 


a 2 sini (Rq) 


n. 


22.87 ± 0.17 


mi sin"* i (Mq ) 


n. 


15.82 ± 0.26 


m 2 sin 3 i (Mq) 


n. 


12.00 ± 0.19 


r.m.s. (km s ) 


4.1 


6.1 



3 The period has been adjusted simultaneously together with the 
other orbital parameters, adopting the values derived in Sect, l3~2l as 
initial guess. 



der of 20° over about 7 yr. This question is addressed in 
the next section. 

3.4. Apsidal motion 

As already mentioned, several authors have collected a 
number of spectroscopic observations of HD 152218 and 
measured the RVs of its t wo components. The largest sets 
are due to IStruvd f)l944h . IHCB74L ISLP97I an dlGMOll who 
respectively obtained 5, 15, 22 and 7 spectra. Struvd pub- 
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Table 5 

Best fit values for the longitude of the periastron u) at different 
epochs. The data sets labelled 'A' to 'C are described in the text. 



0.5 1 

Fig. 7. Literature RV measurements plotted over the He I SB2 
RV-curves of Fig. \E\ shifted according to their respective appar- 
ent systemic velocities. The phases were computed using the He I 
ephemeris of Table [4] Different symbol s indica te different works : 
(five branch star) IStruvel lll944h; (x) lHCB7l (as terisk) ILM83I: 
(open fi ve-branch starl lConti et alj [U977T ); (+l lSLP9 Jj; (Y) corrected 
GMOlJ. These symbols inscribed in open circles indicate secondary 
RV measurements. 

lished the o bserving dates using universal times with day 
and month. IHCB74I later presented corresponding Julian 
dates (JD). We however note that, though quoted other- 
wise, these dates are not heliocentric Julian dates (HJD) 
but simply JD. We thus recomputed the appropriate HJD 
before using the Struvd data set. A word of caution is also 
requir ed for th e lGMOll observations. As already mentioned 
above. iGMOll quoted erroneous HJDs for part of their ob- 
servationO. We thus corre cted th e corresponding values 
prior to the inclusion of the lGMOll observations. We com- 
pleted these four s ets by an additional measur ement from 
Conti et all (|l977t) and another one by ILM83I . We finally 



added the RVs from the present work (obtained as the aver- 
age of the 7-corrected individual lines measurements) and 
we adopted the systemic velocities derived in Sect. 13.31 

In a first approach, Fig. [7] shows the literature measure- 
ments plotted over the He I orbital solution of Tabic |4j The 
secondary measurements present a larger scatter and we 
will thus focus on the primary data for the rest of this 
section. In Fig. O we also observe systematic shifts be- 
tween the literature measurements and our newly derived 
orbital solution. As suggested in the previous section, this 
most probably comes from the presence of an apsidal mo- 
tion. Fixing the period and the eccentricity at the values 
quoted in Table 0] for the He i-solution, we adjusted the 
different data sets of Tableland derived, for each epoch, 



4 In iGMOll 's Table 2, the HJDs around 2 449 910 are probably 
correct, those around 2 450 590 are overestimated by exactly one day 
for observations obtained after midnight UT. We could not check the 
more recent observations obtained bv lGMOl] at HJDs~2 451 360. 



Data set HJD 


- 2 400 000 u 


(°) r.m.s. 


(kms 1 ) 


Struve (1944) 


31223.5 - 


100 


5.4 


HCB74 


40818.9 


34 


11.7 


SLP97 a 


49627.5 


82 


7.6 


This work 


51652.1 


102 


6.6 


A 


49333.8 


72 


19.9 


B 


50861.8 


95 


11.8 


C 


52621.3 


106 


1.2 



a. The first point of the SLP97" set is separated by about 3300 days 
from the other points of the set and was not included in the fit. 



the best-fit longitude of periastron. These values are plot- 
ted in Fig. El We no t e tha t, using the RV sets of lHCB74 
and Stickland et al.1 (1991), we re-derived the longitude 
values respectively quoted in their proposed orbital solu- 
ti ons. We also bu ilt three subsets of data using the data 
of lGM0ll . lSLP97l and ours. These sets, labelled A, B and 
C in Table El correspond respectively to data obtained be- 
fore HJD = 2 450 000, in-between HJD = 2 450 000 and 
HJD = 2 451 800 and after HJD = 2 452 000. These sets 
are however not independent from the sets labelled 'SLP97' 
and 'This work' in Table [5] as they are formed by the same 
observations. We consider our approach to provide a kind 
of moving average over the time. It can also be seen as a way 
to give more weight to the more recent data, th at prob ably 
provide a better accuracy than the older sets of Struvd and 
IHCB74 . Using a linear regression, we then computed the 
best-fit linear relation that reproduces the observed values, 
the slope of which provides the rate of the apsidal motion. 
We obtained: uj = 9?1 ± 0?8 x 10~ 3 d _1 , corresponding to 
about u> = 3?3 ± 0?3 yr _1 . It is worth to note here, that, 
because of the similar K\ and K2, the shape of the RV- 
curve is litt le altered by a change of uj of tt. B y chance , thus, 
neither the Struve (|l944f ) data nor those of lHGB74l could 
strongly bias the orbital period determination of Sect. 13.21 
Indeed similar results are obtained while including or re- 
jecting these two sets from the period search analysis. 

A word of caution is however necessary concerning the 
rate of the apsidal motion. Indeed, we have also performed 
a similar analysis using our data only. These were grouped 
into five different subsets according to the different observ- 
ing epochs. The apsidal motion is again clearly detected, 
even considering the comparatively smaller time span of 
our spectroscopic campaign. The obtained rate is however 
somewhat slower: uj — 1?4 ± 0?6 yr^ 1 . Accounting in ad- 



dition for the ISLP97I data, we rather obtain uj = 3?0 ± 
0?6 yr _1 . Although there is little doubt left about the ex- 
istence of an apsidal motion, the exact rate at which the 
latter occurs might still be uncertain. This probably reflect 
the difficulty to combine inhomogeneous data sets. 



n 



4. Physical parameters 

4.1. Spectral types and luminosity classes 

The dominant spectral signature of the primary com- 
ponent in the spectrum of HD 152218 is characteristic of 
a lat e O-type star. We adopted the classification criteria 
from IConti ( 1973b ) as adapted to late O-stars bv IMathvsl 



100 - 



(1988). They are based on the equivalent width (EW) ratio 



of the He 1 14471 and He n 14542 lines. Doing this, we only 
considered the EWs measured on the spectra where the two 
components are disentangled. For the primary, we obtained 
a mean \ogW[ = logWA447i — logWA4542 = 0.33 ± 0.06 
which corresponds to a spectral type 09, with spectral type 
08.5 within la. For the secondary, we obtained logW^ = 
0.67 ±0.08. This indicates an 09.7 star, with the 09.5 type 
at l-er. 

To determin e the luminosity c la sses, we adopted the 
criterion from IConti fc Alschuleri (|197ll ) based on the 
EW ratio of the Si IV 14089 and He 1 14144 lines. We re- 
spectively obtained log W" = log Ha4089 — log Wa4144 = 
0.25 ± 0.05 and logW£ = -0.12 ± 0.13. These values 
lead to giant and main sequence luminosity classes for 
the primar y and se c ondar y components respectively. Ac- 
cording to iMathvd (|1988l ). we also measured \ogW{" = 
log WA4388+log Wa4686 = 5. 03±0. 04. This corresponds to a 
supergiant if the primary brightness ratio l\ = jj^- > 0.76, 
to a giant if 0.76 > l\ > 0.43 and to a main sequence star 
otherwise. For the secondary, we get log W^" = 4.22 ± 0.15. 
The respective critical values for li = -k 2 - are this time 

C -ktot 

0.12 and 0.07. Such a low con tribution of the secondary to 
the total flux is improbable. Mathvsl criterion thus indi- 
cates a supergia nt secondary, which is clea rly at odds with 
the resul t of the Conti fc Alschulei ( 197ll ) criterion. How- 
ever, the IMathvsl criterion can be biased if the He II14686 
line is partly filled by emission produced e.g. in a wind 
interaction region, which is possibly the case in HD 152218 
(see Sect. [5]) The comparison of our spectra with the atlas 
of standards (Wal born fc FitzpatricklTl990h seems to favor 



a giant class for the primary star, but it is very difficult to 
reject the main sequence option. It also indicates that the 
secondary is most probably a main sequence star. 



4.2. Optical brightness ratio and evolutionary status 

We roughly estimated the optical brightness ratio based 
on the dilution of the primary and secondary lines in the 
spectrum of HD 152218. For this purpose, we compared the 
mean EWs of the primary and second ary lines with typi- 
cal ("average d ) EWs of respec t ively 09 ( Conti fc Alschuleil 
19711 IContil. Il973at IMathvsl Il988l [1989^ and Q9.7 stars 



(|MathvsL I1988T 1989). We first compared the primary line 
strength with typical 09III EWs. Based on the He 1 114026, 
4144, 4388, 4471, He n 114542, 4686 and Si iv 14088 lines, we 
obtained an averaged brightness ratio of 4 = 0.76 ± 0.05. 
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JD (HJD-2 400 000) 

Fig. 8. Values of the argument of periastron u> computed for different 
epochs. The dotted line gives the best-fit relation, corresponding to 
an apsidal motion of 9?1±0?8 X 10~ 3 d~\ or about 3?3±0?3 yr" 1 . 
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Fig. 9. HD 152218 primary (filled symbols) and secondary (open 
symbol) location in the H-R diagram. Triangle (resp. square) symbols 
indicate an adopted gian t (resp. main se q uence ) luminosity class. The 
evolutionary tracks from lSchaller et al.1 ljl992T l have been overplotted 
(plain lines) together with isochrones (dotted lines) computed for 
ages ranging from 2 to 10 Myr with a step of 2 Myr. 

Using typical 09 V EWs yields a lower valued = 0.68±0.11. 
We however note that the Si iv 14088 line is particularly 
strong in the primary spectrum and is the main cause for 
the large scatter obtained. Rejecting the Si IV 14088 line 
from the comparison with 09V typical EWs, we obtained 
h =0.63 ±0.04. 

Turning to the secondary component, we also com- 
pared the secondary line stren gths to the typical EWs 
reported bv lMathvd (Il988l . ll989h . who only considered the 
HeiAA4388, 4471 and HeiiAA4542, 4686 lines. We note 
that the He n 14686 line provides results that clearly differ 
from the three other lines, but this is not surprising since 
the He II 14686 line is particularly faint in the secondary 
spectrum. Mean values for the three other lines indicate 
l 2 = 0.34 ± 0.02 and 0.41 ± 0.05 for dwarf and giant star 
respectively. The obtained values for the primary and sec- 



ondary contribution to the observed fluxes are thus in ac- 
ceptable agre ement with each ot her. They also agree with 
the results of lPennv et al. who quoted a primary to 

secondary UV flux ratio of 2.0. In the following, we will 
adopt a primary to secondary brightness ratio correspond- 
ing to = 0.67 ± 0.10, which reasonably well reproduces 
the results of the above discussion. 

4.3. Effective temperatures and radii 

Sung (2005, private communication) reported V = 7.562 
and B — V = .191 for the system. Using R = 3.3 
( Sung et all. Il998l ) , adop ting the absolute colors quoted 
bv ISchmidt-Kaler ( 1982 ) and a distance modulus DM = 
11.07 ± 0.04 (|Sana et al. . 2006b . and references therein), 
preliminary computations indicate that the system is too 
faint to harbour a giant primary, unless its contribution l\ 
to the observed flux brightness were very close to unity. 
This is clearly not confirmed by the previous discussion and 
we suggest that the primary in HD 152218 could be a class 
IV object. Correcting from the reddening and the distance 
effects, we obtained My,i = —4.72 ± 0.22 and My, 2 = 
—3.95 ± 0.37. Because the intrinsic colour (B — V)q of an 
O-type star does almost not depend on the stellar gravity, 
we note that these values are free from the hypotheses on 
the lu minosity class. According to lHumphrevs fc McElrov 
( 19841) . these values correspond to the primary being a 
V-IV star and the sec ondary, a V star. Adopting the 
Humphreys fc McElrovl temperature scale and bolometric 
corrections, we derived Mboi,i = —8.0 ±0.2 (resp. Mboi.i = 
—7.8 ± 0.2) if the primary is a dwarf (resp. giant) star. We 
obtained Mboi,2 = —7.1 ± 0.4 for the main sequence sec- 
ondary. Finally, using temperatures of T° S 9IU = 32.8^JkK 
(resp. T£g V — 34.0g2;8 kK), we derived a primary radius 
of 10.3 ± 1.3 R (re sp. 9.9 ± 1.2 R w ). For t he secondary 
we interpolated the I Humphreys fc McElrov temperature 



scale to obtain values of T e % 9 - 7V = 31.2ig;g kK. This corre- 
sponds to i?2 = 7.9 ± 1.7 R Q , where the quoted confidence 
intervals have been computed by adopting the parameters 
of the neighbouring spectral sub-types. 

4.4. Masses and orbital inclination 

Fig. [9] presents the location of the stars in the H-R 
diagram. The evolutionary age of both components is 
about 4 to 5 Myr, which is in agreement with previ- 
ous cluster age de t erminations from photometric studies 
dPerry et al. . 1990t Balona fc Lanev . 1995 ; Raboud et al 



1997t iBaume et all I1999F From a rough interpolation of 



the tracks, the main sequence primary had an initial mass 
of about 25 M , corresponding to a present mass close to 
24.5 M Q , slightly larger than typical 09V masses (which 
suggests again that the primary component has reached an 
intermediate stage of evolution between classes V and III) . 
Similarly, the secondary initial mass was about 18.5 M . 
Its present mass should be about 18.2 M Q . 



Comparing the obtained minimal mass of Table 2] with 
typical 09V masses suggests an orbital inclination i of 
about 70°. Comparing with 09III masses would rather 
yield i ~ 60°. Using the constraints deduced on the com- 
ponents radii, we found that the limiting inclination for 
which eclipses may or not occur in the system lies around 
55-65°. Clearly, HD 152218 is a limiting case. Photometric 
monitoring is required to search for the signature of eclipses 
or ellipsoidal deformation that could help to constrain the 
inclination of t he system. Finally, we used the formula of 
lEggletonl (|l983h to estimate the Roche lobe radii. We found 
that, according to the previous results, both stars should 
presently fill about 17% of their Roche lobe volume. 

5. XMM-Newton observations 



5.1. Light curves and spectra 

With a typical count rate of about 12 x 10~ 3 cnts -1 in 
the two MOS instruments, HD 152218 is a relatively bright 
X-ray source. Because the count rates in the hard energy 
band (2.5-10.0keV) are very low, thus associated with large 
uncertainties, we restrain our analysis to energies below 
2.5 keV. Fig. [10] presents the obtained broad band light 
curves in the total band and in the Sx [0.5-1.0 keV] and Mx 
[1.0-2.5 keV] bands. Except for Obs. 2 at (f> w 0.7, the two 
instruments give very consistent results, indicating a slight 
modulation with an increase of the X-ray flux of about 30% 
at apastron compared to periastron. 

To constrain the physical properties of the emitting 
plasma, we adjusted a series of optically thin thermal 
plasm a (mekal) models using the XSPEC software (jArnaudl . 
1996]). We simultaneously fitted the two MOS spectra of 
each pointing, adopting an interstellar column of absorbing 
matter (neutral hydrogen) of -/Vh.ism = 0.29 x 10 22 cm~ 2 . 
This value wa s obtained using E (B — V) = 0.50 and 
the formula of iBohlin et al. ( 19781 ). Except for Obs 2 and 
4, for which the effective exposure time was reduced by 
high background events, better results are obtained using 
two-temperature (2-T) mekal models. Table [6] summarizes 
the best-fit values and gives the intrinsic fluxes and lumi- 
nosities. There is no obvious modulation of the spectral 
parameters. The modulation of the hardness ratios (Fig. 
TU|) however suggests that the observed emission is slightly 
harder when the flux level is high, and lower otherwise. 
The merged spectra are well described by 2-T models with 
kTi = 0.31 keV and kT 2 = 0.71 keV. Fig. [TT] displays 
these spectra together with the best-fit 2-T models. It re- 
veals that a hard emission component (> 5 keV) might be 
present, but clearly the SNR is far insufficient to reliably 
constrain its properties. 



5.2. HD 152218 wind properties 

We used the constraints on the physical properties de- 
duced in Sect. 2] together with the mass loss recipes of 



Table 6. Best-fitting models and resulting X-ray fluxes at Earth. The 1-T model used was wabsisM * wabsC/Vn) * mekal(kT) while the 2-T model corresponds to wabsigM * (mekal(kTi) 
+ wabs(JVH.2) * mekal(kT2 ) ) . The term wabsjgM was fixed to the interstellar value (JVjj ism = 0.29 X 10 22 cm -2 ). Columns 1 and 2 give the phase and the observation number. Columns 
3 to 5 provide the best-fit parameters: TVh, the absorbing column; kT, the model temperature; norm, the normalisation factor (norm = ^p—p J n c nndV with d, the distance to the 
source — in cm — , n c and njj, the electron and hydrogen number densities - in cm -3 ). The quoted upper and lower values provide the limits of the 90% confidence interval. Column 6 
lists the corresponding reduced chi-square and the associated number of degrees of freedom (d.o.f). Columns 7 to 10 provide the observed fluxes in the 0.5 - 10.0 keV energy band and 
in the Sx (0.5 - 1.0 keV), Mx (1.0 - 2.5 keV) and Hx (2.5 - 10.0 keV) bands respectively while Columns 11 to 14 provide the same fluxes, but corrected for the adopted ISM absorbing 
column. Column 15 finally provides the X-ray luminosities in the 0.5 - 10.0 keV band, assuming DM = 11.07. 
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Fig. 10. Broad band X-ray light curves in the two EPIC MOS instruments. The different energy ranges considered are indicated on top of 
each panel. The vertical bars give the ±l-cr uncertainty on the background-corrected count rates. The horizontal bars indicate the duration 
of each pointing. The bottom right panel displays the evolution of the source hardness ratio with the phase, as defined on top of the panel. 



Vink et al. ( 2000l . 2001 ) to estimate the properties of the 
winds in HD 152218. We obtained v^i = 2280 kms" 1 , 
Voo 2 = 2340 kms" 1 , log Mi = -6.59 (M© yr^ 1 ) and 
log A/2 = —7.37 (Mq yr _1 ). The estimated terminal ve- 
locities are only slightly larger tha n the 2050 kms -1 value 
measured bv lHowarth et al.l (|l997l ) . Using these parameters 
with a /3 = 1 acceleration law for the winds, we could not es- 
tablish a ram pressure equilibrium at all phases, indicating 
that the wind collision structure could be quite unstable. 
We suspect that the wind veloc ity law might be s trongly 
affected by radiative inhibition ( Stevens fc Pollock , 19941 ). 
The latter may even govern the flux of colliding matter. 
Clearly, only a detailed modeling could allow to get an in- 
sight into the wind collision in this system. This is far be- 
yond the goal of the present paper. We will thus limit our- 
selves to the observational result that the flux is slightly 
variable and that the maximum of emission is observed 
near apastron while the minimum is reached at periastron. 
As HD 152218 is a close binary, the wind interaction will 
inevitably take place within the acceleration region of the 
winds. We thus expect the collision to be stronger at apas- 
tron, because of the higher pre-shock velocity reached by 
the wind at this particular phase. However, we caution that 
more subtle aspects (such as radiative inhibition) probably 
affect the strength of the winds, and might thus govern the 
actual position and shape of the on-going wind interaction. 



Merged spectrum 




nergy (keV) 



Fig. 11. Top panel: Merged EPIC MOS1 (black) and merged EPIC 
MOS2 (light grey - hardcopy edition - or red — electronic edition) 
X-ray spectra of HD 152218 and best-fit 2T models. Lower panel: 
Individual contributions of the different energy bins to the \ 2 . The 
contributions are carried over with the sign of the deviation (in the 
sense data minus model). 



6. Summary 

We have presented the results of a long-term spectro- 
scopic monitoring campaign on the massive binary HD 
152218. The last two investigations on the object were pre- 



senting conflicting results. The present data complete the 
existing observational sets and allow to definitely overcome 
the ambiguity on the orbital period of the system, bringing 
to firm ground the period of 5.604 d. From our analysis, we 
rej ect the shorter period and h igher eccentricity proposed 
bv lGarci'a fe Mermilliodl (l200lh and rather confirm the ear- 
lier work of Stickland et al.l l 1997 ). based on IUE observa- 
tions. Our newly derived orbital solution is characterized by 
a slightly lower eccentricity than previously accepted: e = 
0.259 ± 0.006. The system is most probably formed by an 
09 sub-giant and an 09.7 main sequence star. We derived 
minimal masses of 15.82 ± 0.26 M and 12.00 ± 0.19 M Q 
and we constrained the component radii to values of R± = 
10.1 ± 1.0 R and R 2 = 8.9 ± 2.0 R Q . These values are 
in good agreement with previous findings and further indi- 
cate that HD 152218 should have an orbital inclination of 
about 60-70°. This corresponds to a limiting case in which 
eclipses or ellipsoidal variations might be observed. 

We also report the results of the monitoring of the sys- 
tem in the X-rays thanks to XMM- Newton observations of 
the cluster. The averaged X-ray spectrum is relatively soft. 
It is well reproduced by a 2-T mekal model with compo- 
nent temperatures about 0.3 and 0.7 keV. We showed that 
the system presents an increase of its X-ray flux of about 
30% near apastron compared to periastron. We note that 
this could be the signature of an ongoing wind-wind inter- 
action process occurring within the wind acceleration re- 
gion. Such a scenario indeed predicts a stronger shock, and 
thus a larger X-ray emission from the colliding zone, when 
the separation between the components is larger, allowing 
the winds to reach larger pre-shock velocities. This is also 
supported by the apparent modulation of the hardness ra- 
tio, which seems higher when the emission level is stronger. 
Preliminary computations however indicate that no simple 
modeling, reproducing the main properties of the interac- 
tion, can be obtained. By contrast, the same computations 
rather point towards second order effects, such as radiative 
inhibition, to play a major role in this system. These effects 
could be strong enough to govern the structure and shape 
of the interaction. HD 152218 might thus be an excellent 
case to test and/or refine wind collision models, including 
detailed physics on the wind acceleration and inhibition. 
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